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Abstract 

Relativistic nuclear collisions data on two-particle correlations exhibit structures as function of relative azimuthal angle and rapidity. 
A unified description of these near-side and away-side structures is proposed for low to moderate transverse momentum. It is based 

CD on the combined effect of tubular initial conditions and hydrodynamical expansion. Contrary to expectations, the hydrodynamics 
solution shows that the high energy density tubes (leftover from the initial particle interactions) give rise to particle emission in two 
directions and this is what leads to the various structures. This description is sensitive to some of the initial tube parameters and 
may provide a probe of the strong interaction. This explanation is compared with an alternative one where some triangularity in the 

^> [initial conditions is assumed. A possible experimental test is suggested. 

CN Keywords: Relativistic heavy-ion collisions, Particle correlations and fluctuations, Collective flow 
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1. The need for a unified description 

One of the most striking results in relativistic heavy-ion colli- 
sions at RHIC and the LHC, is the existence of structures in the 
two-particle correlations Minna! plotted as func- 
tion of the pseudorapidity difference Arj and the angular spacing 
A<p. The so-called ridge has a narrow A(f> located around zero 
and a long A77 extent. The other structure located opposite has a 
single or double hump in A<f>. In order that two particles, emit- 
ted at some proper time r /.„„,, appear as correlated over several 
rapidity units, the process that correlated them must have oc- 
curred llgl Qii at a much smaller proper time due to causality. 
Therefore, the existence of long range pseudorapidity correla- 
tions must be related to early times in the nuclear collisions and 
thus has motivated many theoretical investigations. 

Hydrodynamics has now been established as a good tool to 
describe many data from relativistic heavy-ion collisions so it 
should be able to provide a description for the above mentioned 
structures (for low to intermediate transverse momenta). In fact, 
as noted with RHIC data, a hydrodynamics based explanation 
is attractive because of the various similarities (see e.g. [11]) 
between bulk matter and ridge (transverse momentum spectra, 
baryon/meson ratio, etc). In addition, it was shown (particularly 
at the LHC) that particle correlations can be understood in term 
of anisotropic flow Fourier coefficients Il2l.ll3l . ll4n . This points 
towards the necessity to have a unified hydrodynamic descrip- 
tion of near and away-side structures. 

In early models, it was suggested that the combined effect of 
longitudinal high energy density tubes (leftover from initial par- 
ticle collisions) and transverse expansion was responsible for 
the ridge [15, 16, 9l[l0lll7|]. The particle emission associated 
to the tube was expected to occur in a single direction, so this 
would cause a ridge but no away-side structure. In addition, the 



effect of hydrodynamics was usually assumed to be of a certain 
type (e.g. a blast wave in 

SElntl) and in fact when hydrody- 
namic expansion was actually computed, it seemed to lead to a 
disappearance of the initial high energy density tubes H 1 811 and 
therefore of their particle emission. 

In a previous work 01911 . we presented evidence that hydro- 
dynamics might in fact reproduce all structures using the NeX- 
SPheRIO code. This code starts with initial conditions from the 
event generator NeXus I20I1 and solves the hydrodynamic equa- 
tions on an event-by-event basis with the method of Smoothed 
Particle Hydrodynamics. In I0, the NeXSPheRIO events 
were analyzed in a similar way to the experimental ones, in par- 
ticular the ZYAM method was used to remove effects of elliptic 
flow. We later developed a different method to remove elliptic 
flow from our data and checked that all structures were indeed 
exhibited and other features well reproduced (dependence on 
the trigger- or associated-particle transverse momentum, cen- 
trality, in-plane/out-of-plane trigger, appearance of a peak on 
the ridge). However, when using NeXSPheRIO, it is not clear 
how the various structures in the two particle correlations are 
generated. The aim of this paper is to investigate this, study- 
ing in detail what happens in the vicinity of an energetic tube 
(section 2) and then extending the results to a more realistic 
complex case (section 3). We will also compare our explana- 
tion with an alternative one that assumes some triangularity of 
the initial conditions (section 4). 

2. A simplified model 

2.1. Origin of the near-side and away-side structures 

We will consider central collisions only and use a simplified 
model. Fig. [T](left and center) shows a typical example of ini- 
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Figure 1 : Left and center: initial energy density for a NeXus central Au-Au collision at 200 GeV A in the y = and r\ = plane respectively. Right: Comparison of 
the parametrization given by Eq.(l) (solid lines) with the original NeXus energy density (dashed lines), along the lines a-e (in the r\ = plane). 



tial conditions (initial energy density) obtained in NeXus with 
various tubular structures along the collision axis. 

The origin of these structures is the following. To model soft 
physics in p-p collision, it is common to assume that strings (or 
color flux tubes) are formed, either via the excitation of the pro- 
tons or due to color exchange between them. In A-A collisions, 
these strings may overlap leading to longitudinally extended re- 
gions of higher energy density such as those in Fig. [1] An alter- 
native description of A-A collisions, based on gluon saturation, 
is that the two colliding nuclei can be viewed as Color Glass 
Condensates. Shortly after their collision, these produce strong 
color flux tubes called "Glasma". Therefore the possibility that 
tubular structures exist in the initial conditions is general but 
their exact characteristics are not known. 

In the simplified model, only one of the high energy tubes 
from NeXus (chosen close to the border) is considered and the 
complex background is smoothed out. This leads to the follow- 
ing parametrization of the initial energy density 
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where ro = 5.4 fm. A comparison of this parametrization with 
the original NeXus energy density is shown in Fig. Q] (right). 
Except for the inner region (which has little importance cf. 
§2.2), the agreement is reasonable. We use this parametriza- 
tion in order to have a realistic tube description. However as 
already mentioned, the exact characteristics of the color tubes 
are not well known, therefore later we will consider various 
variations of the parameters. 

In this simplified model or one-tube model, transverse ex- 
pansion is computed numerically while longitudinal expansion 
is assumed boost-invariant, until freeze out at some constant 
temperature. The resulting single-particle angular distribution, 
shown in Fig. |2](top), has two peaks located on both sides of the 
position of the tube (placed at = 0) with separation ~ 2 (this 
is not a parameter), more or less independently of the value of 
the transverse momentum. Particle emission is computed as- 
suming sudden freeze out. Since this is an approximation to 
real particle emission, we have checked that varying the freeze 
out temperature (between 130 and 150 MeV) does not affect 
qualitatively our result. 
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Figure 2: Angular distributions of (direct) charged particles in some different 
pr intervals (top) and resulting two-particle correlations (bottom) in the sim- 
plified model (for a freeze out temperature of 0.14 GeV). 



This two-peak emission is in contrast with what happens 
when a blast wave is assumed, namely the fact that high-energy 
tubes emit in a single direction. However, its occurrence can be 
understood from Fig. [3] As time goes on, as a consequence of 
the tube expansion, a hole appears at the location of the high- 
energy tube (as in fll8IO . This hole is surrounded by matter that 
piles up in a roughly semi-circular cliff of high-energy-density 
matter, guiding the flow of the background matter into two well- 
defined directions. The two extremities of the cliff emit more 
particles than the background, this gives rise to the two-peaks 
in the single-particle angular distribution. The emission is not 
quite radial as shown by Fig. [3](right), indicating that there was 
a deflection of the background flow due to the pressure exerted 
by the high-energy tube. As seen in fig. 3, the fluid velocity is 
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Figure 3: Temporal evolution of energy density for the simplified model at times 1.0, 3.5 and 8.5 fm (thicker outer black curve corresponding to freeze out 
temperature of 0.14 GeV). Arrows indicate fluid velocity on the freeze out surface (vector length equivalent to 2 fm corresponds to light speed). 



larger at the two extremities of the cliff and smaller nearby, this 
is why in fig. 2 the angular distribution is narrower for larger p, 
particles. 

From Fig. |2](top), we can guess how the two-particle an- 
gular correlation will be. The trigger particle is more likely 
to be in one of the two peaks. We first choose the left-hand 
side peak. The associated particle is more likely to be also in 
this peak i.e. with Acp — or in the right-hand side peak with 
A(p ~ +2. If we choose the trigger particle in the right-hand 
side peak, the associated particle is more likely to be also in 
this peak i.e. with A<p = or in the left-hand side peak with 
A<p ~ -2. So the final two particle angular correlation must 
have a large central peak at A<p = and two smaller peaks re- 
spectively at A<p ~ +2. Fig. |2](bottom) shows that this is indeed 
the case. The peak at A0 = corresponds to the near-side ridge 
and the peaks at A<p ~ +2 form the double-hump ridge. We 
have checked [21] that this structure is robust by studying the 
effect of the height ( 1 2 + 3 in the first term on the right hand side 
of eq. (1)) and shape of the background (r 5 to r 2 in the same 
term), overall initial transverse velocity (increasing radially up 
to 0.6), height, radius and location of the tube (some details are 
shown in the next section). The model was also generalized 
to non-central collisions and the in-plane/out-of-plane trigger 
dependence studied [21]. 

2.2. Dependence on the tube parameters 

In the above calculation, the tube extracted from NeXus ini- 
tial conditions has a radius Ar of order 0.9 fm and (maximum) 
energy density e, of order 12 GeV fm~ 3 (at proper time 1 fm), 
as can be seen from the second term on the right hand side of 
eq. (1) or from fig. 1. Changing Ar affects the height of the 
peaks and spacing in the two-particle correlation as shown in 
Fig. S](top). Changing e, has a similar strong effect. On the 
other side, if Ar and e, are changed maintaining constant the 
energy per unit length E, oc e, (Ar) 2 , as shown in Fig. [4] (bot- 
tom), the two-particle correlation maintains its overall shape, 
the angle between the peaks is almost unchanged and the peak 
heights change less. Therefore a good parameter to character- 
ize the two-particle correlation is the tube total energy per unit 
leng th and tubes thinner than 0.9 fm are not excluded (see also 
[22J]). For this comparison the background was kept unchanged, 



so the thinner tube energy density is really much higher than the 
background one, more realistic cases could be studied but this 
goes beyond the scope of this paper. 
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Figure 4: Top: two-particle correlation for tubes with different radius but sim- 
ilar energy density. Bottom: two-particle correlation for tubes with different 
radius but similar energy content. (Tube position fixed at ro = 5.4 fm.) 

The observation that the correlation is characterized by the 
energy per unit length E, (not by Ar or e, separately) is con- 
sistent with the fact that what matters is the amount of energy 
available for the tube to push the surrounding matter. This also 
explains why the precise shape of the tube energy density (e.g. 
gaussian as in eq. (1)) is not crucial. 

Finally, in Fig. [5] the tube is located (above the background) 
at various distances ro from the center, the height of the two- 
particle correlations at A<p = saturates between 5.5 and 8 fm 
and decreases strongly for smaller distances to the center. 

The physical reason for the behavior of the maximum height 
of the two-particle correlation as function of the tube position 
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Figure 5: Height of the two-particle correlation at 
tube position with respect to the center. 



as function of the 



can be understood by looking at the temporal evolution of the 
energy density. When the tube is exactly at the center, there is 
no privileged direction of emission. When the tube is close to 
the outer border two privileged directions of emission appear 
(cf. Fig. [3]). When the tube is at some small distance such as 
2 fm, even though the strong expansion of the tube presses the 
surrounding matter creating a hole, this happens too much in- 
side to cause the appearance of the two privileged directions of 
emission. We conclude, as mentioned earlier, that only periph- 
eral tubes are important for the particle correlatiorQ. 



3. More realistic case 

With these information (and using the same two-dimensional 
hydrodynamic model as in previous section), we can discuss 
what happens in a more complex event such as a NeXus event. 
In such an event, only the outer tubes are expected to be rel- 
evant, for example in fig. [6] (left), we can pinpoint five such 
tubes, indicated by crosses. When the time evolution of this 
matter slice is studied, holes appear in the vicinity of the for- 
mer location of the tubes indicated by crosses in fig. |6](center). 
Due to expansion and the fact that one tube can interfere with 
another, we do not expect perfect one-to-one correspondence 
(though in this particular event, it is approximately the case). 

The shape of the two-particle correlations for a single tube 
(in particular the peak spacing) is relatively independent of its 
features so the various tubes in the NeXus event under study 
will contribute with rather similar two-peak emission pattern at 
various angles in the single-particle angular distribution. As a 
consequence, the two-particle correlation of this NeXus event is 
expected to have a well-defined main structure similar to that of 
the one-tube model of the previous section (Fig. [2]i surrounded 
by several other peaks and depressions due to trigger and as- 
sociated particles coming from different tubes. This is indeed 
the case as shown in Fig. [6] (right) . The additional peaks and 
depressions have positions depending on the angle of the tubes 



An opposite conclusion was reached in fig. 10 of [23] . A possible expla- 
nation for this discrepancy, is that they use the Gubser solution which has some 
unphysical features at large radii as discussed in 1 24] (in particular see fig. 3 
there). 



between them. When averaging over many events, these inter- 
ference terms cancel out and only the main one-tube-like struc- 
ture is left [19]. In other words: the picture derived in section 2 
also applies to more complex events such as Nexus ones. 



4. Conclusion and perspectives 

Usually, the initial conditions in the hydrodynamic descrip- 
tion of relativistic nuclear collisions are assumed to be smooth. 
It seems however that each time more, understanding data re- 
quires a knowledge of the fluctuating event-by-event initial con- 
ditions rather than an assessment of some adequate smooth ini- 
tial conditions: fluctuations in elliptic flow II25I1 (perhaps the 
very behavior of elliptic flow as function of pseudorapidity 
||26[| ). Fourier coefficients of the azimuthal anisotropic flow (see 
below), two-particle correlations (see introduction). 

In this paper, a unified picture for the structures observed in 
two particle correlations at low to moderate transverse momen- 
tum has been presented. It is based on the presence of longitu- 
dinal high energy density tubes in the initial conditions. These 
tubes are leftover from the initial particle interactions. During 
the hydrodynamic evolution of the fluid, the strong expansion of 
the tubes located close to the border piles up matter in two sym- 
metrical directions, leading to two-particle correlations with a 
near-side and a double hump away-side ridges (for central col- 
lisions). 

An alternative unified picture has also been suggested, the 
idea is the following Ir27ll . The ellipticity of the interaction re- 
gion in a collision gives rise to elliptic flow because of the larger 
pressure gradients along the minor axis of the ellipse. Similarly, 
if the interaction region has some triangular shape, this causes 
triangular flow (due to existing larger pressure gradients in cer- 
tain directions). Both near-side and double hump away-side 
ridges are a natural consequence of triangularity and triangular 
flow 10. 

More generally, it has been suggested that the initial trans- 
verse density in the overlap region can be decomposed using 
an infinite set of moments and any observable (in particular 
anisotropic flow Fourier coefficients v„) can be written as a 
function of these moments [ 2811 . Conventional eccentricities 
such as ellipticity and triangularity, are basically a subset of 
these moments, which may or not be sufficient to characterize 
the v„'s (see e.g. Ir29ln . 

Both tube configuration such as in this paper and triangular- 
ity ll27ll (or more general geometrical shape of the interaction 
region) lead to non-zero eccentricities. In the case of NeXus, 
these eccentricities reflect the initial conditions where a lot of 
peaks or "hot spots" and valleys are present as can be seen in 
Fig. 1 and 6 (similar features can be seen in EPOS initial condi- 
tions i22ll ). On the other side, depending on the sampling pro- 
cess for the nucleon-nucleon collisions and sources of fluctua- 
tions included, the interaction region may present a geometrical 
shape without dominance of a few tubes and might be described 
by conventional eccentricities such as suggested in 12711 . 

The hydrodynamic al response to the anisotropics differs in 
the two cases. In our description, the structures in the two- 
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Figure 6: Left: initial energy density for a NeXus central Au-Au collision at 200 GeV A in the Tj = plane with five outer tubes. Center: appearance of holes in the 
energy density during the temporal evolution (obtained assuming longitudinal boost invariance). Right: two-particle correlations. 



particle correlations and the various v„'s are a response to in- 
dividual outer tubes: this is a local effect. In the case of trian- 
gularity (or more general geometrical shape) of the interaction 
region, the structures in the two-particle correlations and v„'s 
correspond to the various geometrical deformations: this is a 
global effect. In our approach, the angular size of the near- 
side ridge, ~ 2, is of local hydrodynamic origin and can be 
roughly understood as the diameter of the hole divided by the 
radius of the nucleus, i.e. 2c s Tf. 0Ut /R (ignoring the non-radial 
flow of particle emission cf. §2.1). In the triangularity case, 
this angular size is fixed by global geometry (for example 2n/3 
in the simple case of an equilateral triangular shape). Simi- 
larly, in our approach, the relative height of the peaks in the 
two-particle correlation is of local hydrodynamical origin: for 
example, for a single tube (in a central collision), the highest 
peak is approximately twice higher than the two smaller peaks 
(cf. fig. 2). In the triangularity case, the relative height of the 
peaks in the two-particle correlation is related to global geom- 
etry: for example, in the simple case of an equilateral trian- 
gular shape, there should be three equal height peaks in the 
two-particle correlation. However, to turn the connection be- 
tween initial (local or global) geometry and flow more precise, 
pre-equilibrium evolution B30I1 and viscosity 113 ill should be in- 
cluded since both, though they are small effects, can smear out 
the initial anisotropies. 

Finally, to further test the presence of tubular structures in 
the initial conditions, we suggest to build 2+1 correlations, fix- 
ing both a trigger particle and a first associated particle, this 
last one with A<f>i ~ 2. This choice ensures that both particles 
come from different emission peaks in the single particle angu- 
lar distribution. Then, in our approach, the second associated 
particle will be more likely to come from the same emission 
peak as the trigger, i.e. with A<p2 — or the same emission peak 
as the first associated particle, i.e. with A<p2 ~ 2. Naively, the 
plot of the 2+1 correlation as function of Afc vs. A772 should 
present two stripes located at A<p2 = and A<p2 ~ 2. How- 



ever, in practice, there appears a third weaker stripe, due to the 
background. In the eccentricity case, the plot is expected to 
be more complicated (in the simple case of only equilateral tri- 
angular shapes - with longitudinal extension, there should be 
three equally bright stripes). Additional work with realistic ini- 
tial conditions is needed to check whether 2+1 correlations may 
indeed permit to distinguish between both scenarios. 

Since the two scenarios, isolated tube configuration such as 
in this paper and triangularity (or more general geometrical 
shape of the interaction region) correspond to different schemes 
of initial state energy deposition, distinguishing between them 
may allow to learn about how the strong interaction proceeds 
during high energy nuclear collisions. 
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